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SUMMARY 


A unidirectional fiberglass epoxy composite specimen is modelled 
as a homogeneous transversely isotropic continuum plate -like medium 
Acousto-ultrasonic non-contact input-output characterization by trac- 
ing the P waves in the continuum is studied theoretically with a 
transmitting and a receiving transducer located on the same face of 
the plate The isotropic plane of the equivalent continuum plate 
model lies in the midplane of the plate and is parallel to the top and 
the bottom faces of the plate. 

The single reflection problem for an incident P wave at a stress - 
free plane boundary in a semi- infinite transversely isotropic medium 
whose isotropic plane is parallel to the plane boundary is analyzed 
first. It is found that an obliquely incident P wave results in a 
reflected P wave and a reflected SV wave. One quadrant of each of two 
sheets of the slowness surfaces of the reflected P and SV waves is 
plotted The angle of reflection of the reflected P wave is equal to 
the angle of incidence of the incident P wave. However, the angle of 
reflection of the reflected SV wave is smaller than the angle of inci- 
dent of the incident P wave. Consequently, no critical angle phenome- 
non occurs. The amplitude ratios of the reflected P and SV waves to 
the incident P wave as a function of the angle of incidence are 
plotted. The balance in energy flux normal to the plane boundary is 
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checked. 


The delay time for propagation between the transmitting and the 
receiving transducers is computed as if the P waves were propagating 
m an infinite half space. It is found that the directional depen- 
dence of the phase velocity of the P waves propagating in the trans- 
versely isotropic medium has a significant effect on the delay time, 
as opposed to the directional independence of the phase velocity of 
the P waves propagating in an isotropic medium. 

The displacements associated with the P waves in the plate and 
which may be detected by the non-contact receiving transducer are 
approximated by an asymptotic solution for an infinite transversely 
isotropic medium subjected to a harmonic point load. The polar dia- 
grams for the directivity functions are plotted at frequencies of 
0 75, 1 50 and 2.25 MHz. 

The study enhances the quantitative understanding of acousto- 
ultrasonic nondestructive evaluation (NDE) parameters such as the 
stress wave factor (SWF) and wave propagation in fiber reinforced com 

posites or any other materials which can be modelled as transversely 
isotropic media 
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INTRODUCTION 


Fiber reinforced composite materials are attractive materials for 
aerospace applications because of their high specific mechanical prop- 
erties. It has been shown that many composites, such as fiberglass 
epoxy composites or fiber reinforced ceramics, as shown in Fig 1, may 
be modelled as homogeneous transversely isotropic continua [1] In 
this work, acousto-ultrasonic (AU) non-contact input-output character- 
ization of a homogeneous transversely isotropic elastic plate is 
investigated by tracing P waves. 

First, the single reflection problem of an incident P wave at a 
stress -free plane boundary in a semi- infinite transversely isotropic 
medium, whose isotropic plane is parallel to the plane boundary is 
considered. At such boundaries, the conditions for the existence of 
wave mode conversion, the angle of reflection of the reflected wave, 
and the amplitude ratios of the reflected wave to the incident wave 
are derived. 

Second, the P wave input-output relations are derived when mul- 
tiple reflections occur at the top and bottom faces of the plate. The 
delay time between input and output versus the distance separating the 
transmitting and receiving transducers is analyzed. The directivity 
functions of the stresses associated with the P waves are computed 
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And, the output displacement at the non-contact receiving transducer 
is approximated by an asymptotic Solution. 

This investigation should enhance the quantitative understanding 
of AU NDE parameter such as the stress wave factor. It also provides 
the potential for assisting in the development of better NDE schemes 
utilizing the SWF as well as other AU parameters. 
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SINGLE REFLECTION PROBLEM AT STRESS -FREE PLANE BOUNDARY IN SEMI- 


INFINITE TRANSVERSELY ISOTROPIC MEDIUM WHOSE ISOTROPIC PLANE IS PARAL- 
LEL TO PLANE BOUNDARY FOR INCIDENT P WAVE 


1 . Reflected P and SV Waves 


For a homogeneous linearly elastic transversely isotropic contin- 
uum, the number of independent elastic constants is five [1]. Define a 
coordinate system (x, y, z) for a semi- inf inite transversely isotropic 
medium whose isotropic plane is parallel to the plane boundary where 
the reflection occurs as follows: the plane boundary contains the x 

and y axes, and the z axis is the zonal axis of the medium, which is 
in the direction parallel to the fiber direction shown in Fig. 1. See 
Fig. 2. The generalized Hooke's law is written, relative to the (x, 
y, z) coordinate system, as [1] 


r xx " c ll u > x + c 12 v *y + c 13 w > z 
r yy “ c 12 u - x + c ll v >y + c 13 w > z 

r zz “ c 13 u * x + c 13 v -y + C 33 w,z (1) 

r xz “ c 44( u > z + w . x ) 
r yz “ c 44( v > z + w .y) 
r xy “ c 66< u >y + v . x ) 

where r rs (r, s ■ x, y and z) are the normal (r - s) and shear (r * s) 
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stresses with respect to the coordinate system (x, y, z) ; u, v and w 
are the displacement components of a point in the medium along the x, 
y and z axes, respectively; denotes partial differentiation with 
respect to the variable which follows; and C^, C^, 0^3, C33 and C44 
are the five independent elastic constants where Cgg - 1/2 (C^i - C^) 


Let a plane progressive wave be represented as [2] 

(u, v, w) - A (P x , P y , P z ) exp{iw(S x x + S y y + S z z - t)} (2) 

where S x , S y and S z are the components of the slowness vector, which 
is in the same direction as the normal to the wavefront and whose mag 
nitude is equal to the reciprocal of the magnitude of the phase velo- 
city [1], along the x,y and z axes, respectively, P x , P y and P z are 
the components of a unit vector of particle displacement along the x, 
y and z axes, respectively; A is the amplitude of particle displace- 
ment; t denotes time and w denotes radian frequency. 


It follows from Eqs. (1) and (2) that the stresses 
sented as 

r xx - iwA[CnS x P x + C^2^ y Py + C^3S Z P Z ] exp { iw(S x x + S y y + 
Tyy - iwA[Ci 2 S x P x + C^SyPy + C 1 3 S z P z ] exp { i w ( S x x + S y y + 
r zz - iwA[C 13 S x P x + C 13 S y Py + C 33 S Z P Z ] exp{ iw(S x x + S y y + 
t xz - iwA[C44S z P x + C44S x P z ]exp{iw(S x x + S y y + S z z - t) } 
r yz - iwA[C44S z P y + C44S y P z ] exp{ iw(S x x + S y y + S z z - t) } 
r xy - i- wA [ c 66 s y p x + C 66 S x P y ] exp{ iw(S x x + S y y + S z z - t) ) 


can be repre 

S z z - t)} 

S z z - t)} 

S z z - t)} 

(3) 
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The stress boundary conditions on the stress -free plane boundary 


require that [2] 


(I) 

+ 

r (R) 

xz 


r xz 

(I) 

+ 

r (R) 

yz 


r yz 

(I) 

zz 

+ 

r (R) 
T zz 


(4) 


where r rz ^) (r-x, y and z) represents stresses on the plane boundary 
associated with the incident P wave, and r rz ^) (r - x, y and z) 
represents stresses on the plane boundary associated with the 
reflected waves. 


In order to satisfy Eq. (4), 
w, of the reflected waves be equal 
that [ 2 ] 


it is required that the frequency, 
to that of the incident wave and 


S X (I) - s x W 

(5) 

S (I) - S < R ) 

kJy kJy • 

As a result of Eqs. (5), the slowness vectors of the incident 
and reflected waves lie in a plane called the plane of incidence. 

This analysis can be simplified by assuming that the plane of inci- 
dence is the x - 0 plane; that is, the slowness vectors of the inci- 
dent and reflected waves are in the x - 0 plane, as shown in Fig 2 
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Then it follows from Eqs. (5) that 


c (I) _ q (R) _ A 

^x" “ s x ' “ °- (6) 

It has been shown [3] that P waves or SV waves travelling in a 
plane containing the zonal axis, z axis, of a transversely isotropic 
medium are quasi-longitudinal and quasi- transverse , respectively; that 
is, the components of the unit vector of particle displacement for 
either the P or SV waves along both the y and z axes, P v and P do 

y 2 

not vanish; whereas the components along the x axis, P x , do vanish 
Therefore, it follows from Eqs. (3) and (6) that stresses associated 
with P and SV waves are 

r xy " r xz - 0 

(7) 

r yz * 0 >• T xx * 0 : T yy * 0 ; r zz * 0 . 

It has also been shown that an SH wave travelling in a plane 
containing the zonal axis, z axis, of a transversely isotropic medium 
possesses a transverse displacement only; that is, for the coordinates 
i- n Fig- 2 (P x , Py, P z ) — (1,0,0) [3] Therefore, it follows from 
Eqs. (3) and (6) that the stresses associated with SH waves are 


r xx “ T yy “ T zz “ T yz “ ® 
r xz ** 0 '< r xy * 0. 
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Assume that a P wave is incident on the plane boundary It fol- 


lows from Eqs. (4) and (7) that 


'yz (I) - 0; r zz (D h 0 


' xz 


(I) 


0 . 


(9) 


As a result of Eq.(9), it is immediately known from Eq (4) 
that r xz ^ is equal to zero. This means that no SH wave will be 
reflected back into the medium because a reflected wave of the SH 
type would result in nonzero values of the stress r xz . So, Eqs. (4) 
reduce to: 


r yz 

r zz 


(I) 

(I) 


+ 

+ 


r (R) 

T yz 

r (R) 

1 zz 


0 

0 . 


( 10 ) 


Since either a reflected P wave or a reflected SV wave results 
in nonzero values of the stresses, r y Z ^^ and r Z2 ^. it is therefore 
concluded from Eqs. (7) and (10) that both a P wave and an SV wave may 
be reflected back into the medium. 


2 . Slowness Surface for P and SV Waves 

The equations of motion relative to the coordinate system 
(x,y,z) 
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are [1] 


XX, X 

+ 

r xy ,y 

+ 

r xz , z 

- pu.tt 


xy,x 

+ 

r yy.y 

+ 

r yz,z 

- pv.tt 

(11) 

xz ,x 

+ 

r yz,y 

+ 

r zz , z 

- pw.tt 



where the body forces are identically zero for the homogeneous solu- 
tion. 


It follows from Eqs (1), (2) and (11) that the following equa- 
tions of motion are obtained: 

t c ll s x 2 + c 66 s y 2 + c 44 s z 2 ‘ H p x + (C]_2 + c 66> s x s y p y 
+ (C 13 + C 44 )S x S y P z = 0 

( c 12 + C 66 )S x SyP x + [C 66 S X 2 + C 1;L Sy 2 + C 44 S Z 2 - p]P y (12) 

+ (C^3 + C 44 )SyS z P z “ 0 

(C]_3 + C 44 )S X SyP X + (C 13 + C 44 )SyS Z Py + [ C 4 ^ ( S^j- 2 + Sy- 2 ) 

+ c 33 s z 2 'P] p z ” 0 

The condition for the existence of the plane wave solution is 
expressed by setting the determinant of the matrix of the coefficients 
of P x , Py and P z in Eq. (12) equal to zero [1] : 
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C c 11 S x 2 + C 66 S y 2 + C 44 S i 2 

(C 12 + C 66 )S x S y 

(C 13 + C 44 )s x s z 

< c 12 + C 66> S x S y 

C c 66 s x 2 + c 11 s y 2 + c 44 s z 2 

(C 13 ♦ c 44 )s y s z 

(c 13 ♦ c 44 )s x s 2 

< c 13 + c 44> s y s z 

cc 44 <s x 2 + s y 2 ) + c 33 s z 2 -p] 


(13) 


By expanding Eq. (13), three sheets of slowness surface are 
obtained The slowness surface for an P wave is given in [3] as 


c ll +c 44 2 2 c 44 +c 33 2 

( ) (S x +S ) + ( ) s z 

2 2 

+ 1/2 {[(C n -C 44 )(S x +s5 ) 

+ (C33 -C 4 4) S zl 2 
2 2 2 

+ 4(S x+ S y )S z [(C n -C 44 )(C 3 3-C 44 ) 


-< c 13 +c 44 ) 2 ]) 1/2 - P- (I 4 ) 

Similarly, the slowness surface for an SV wave is given in [3] 


( c ll +c 44 

2 


)(S X +S y ) + (- 


c 44 +c 33 

2 



- 1/2 ( [ (C 11 'C 44 ) (S x +Sy) + (C 33 -C 44 ) S z ] 2 
+ 4(S^+s2)s 2 [(C n -C 44 )(C 3 3 -C 44 ) 
-((C 13 +C 44 ) 2 ] J 1 / 2 “ P • (15 


n 



In the present study, the slowness vectors of the incident P 
wave and the reflected P and SV waves are confined in the y-z plane, 
as shown in Fig. 2. Take numerical values of elastic constants and 
density given in [ 1 ] for the unidirectional fiberglass epoxy composite 
as follows: - 10.581 x 10^ N/m^, - 4.679 x 10^ N/m^, C 33 = 

40 741 x 10^ N/m^, C 44 - 4.422 x 10^ N/m^ , and p - 1850 kg/m^ . One 
quadrant of the intersection of the slowness surface of a P wave 
travelling in the unidirectional fiberglass epoxy composite with the 
plane x ~ 0 and one quadrant of the intersection of the slowness sur- 
face of an SV wave travelling in the unidirectional fiberglass epoxy 
composite with the plane x - 0 are obtained by substituting the numer- 
ical values of the elastic constants and the density given above into 
Eqs. (14) and (15), respectively, and are shown in Fig. 3. 

3 Angle of Reflection 

It follows from Eqs. (5) and ( 6 ) that the y-component of the 
slowness vector of an incident P wave is equal to the y-component of 
the reflected P wave as well as the y-component of the reflected SV 
wave. Accordingly, the relation between the y-component of the slow- 
ness vector of an incident P wave and the y-components of the slowness 
vectors of the reflected P and SV waves is given as 

Sy (i> - SyW - s y (sv) . b (16) 
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where represents the y-component of the slowness vector of an 

incident P wave; Sy^ P ) represents the y-component of the slowness 
vectors of the reflected P wave; Sy(^) represents the y-component of 
the slowness vector of the reflected SV waves; and b is a common con- 
stant 


It follows from Eqs . (6), (14) and (16) that the relation 
between the z-component of the slowness vector of an incident P wave 
and that of the reflected P wave is 

- -S z ( p ) (17) 

The minus sign is due to the fact that the slowness vector of an 
incident P wave points out of the medium, whereas the slowness vector 
of the reflected P wave points into the medium, as shown in Fig. 2 
Consequently, the value of the z-component of the slowness vector of 
an incident P wave, is negative, whereas that of the reflected 

P wave, S z ( p ), is positive. 

The angle of reflection is defined as the angle between the 
slowness vector of a reflected wave, either type P or SV, and the 
normal to the plane boundary where the reflection occurs. Similarly, 
the angle of Incidence is defined as the angle between the slowness 
vector of an incident P wave and the normal to the plane boundary, as 
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shown in Fig. 2. Therefore, the angle of reflection of a reflected P 


wave 0p is 


e v - tan- 1 (S y < p )/S z < p >) (18) 

and the angle of incidence of an incident P wave is 

6 Z - tan- 1 (S y ( I )/-S z < I >). (19) 

It follows from Eqs. (16), (17), (18) and (19) that the angle 
of incidence of an incident P wave is equal to the angle of reflection 
of the reflected P wave, as shown in Fig. 2. 

However, the angle of reflection of the reflected SV wave is not 
equal to the angle of incidence of the incident P wave. For a given 
value of b in Eq. (16), two values of the z-component of the slowness 
vector, S z , of equal magnitude but opposite sign of an SV wave travel- 
ling in the plane x - 0 in the transversely isotropic medium can be 
obtained from the slowness surface for the SV wave from Eqs (6) and 
(15) . The positive z-component of the slowness vector corresponds to 
the reflected SV wave and is denoted as S z ^^. Similarly, for a 
given value of b, there exists a positive z-component of the slowness 
vector, S z ( p >, corresponding to the reflected P wave. Accordingly, 
for any given value of b in Eq. (16), there exist a positive z-compo- 
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nent of the slowness vector for the reflected P wave, S z ^, and a 
positive z-component of the slowness vector for the reflected SV wave, 
S Z CSV). in fact, the z-components of the slowness vectors of the 
reflected P and SV waves, S z ^) and S Z ^V) ( for an incident P wave 
which determines the value of b in Eq. (16) and which travels in the 
plane x - 0 in the transversely isotropic medium, can be obtained from 
the lengths of the perpendicular lines between the abscissa represent- 
ing the value of b and the interesections with two sheets of slowness 
surface for the reflected P wave and the reflected SV wave, respec- 
tively, as shown in Fig. 3. It is apparent from Fig 3 that the 
z-component of the slowness vector of the reflected P wave, S z ^) , is 
less than that of the reflected SV wave, S z ^ sv ^ , for an incident P 
wave travelling in the plane x - 0 in the unidirectional fiberglass 
epoxy composite under consideration [1]. Consequently, the angle of 
reflection of the reflected SV wave 0 S y defined similarly to Eq (18) 
as 

6 SV - tan' 1 (S y < SV >/S z <SV) }> (20) 

is smaller than the angle of reflection of the reflected P wave, from 
Eqs. (16), (18) and (20). The angle of reflection of the reflected SV 
wave as a function of the angle of incidence of an incident P wave 
travelling in the unidirectional fiberglass epoxy composite [1] is 
shown in Fig. 4. 
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The critical angle is defined as the angle of incidence of an 
incident P wave at which the slowness vector of the reflected SV wave 
(or P wave) becomes tangent to the plane boundary where the reflection 
occurs For the case of an incident P wave travelling in the unidirec- 
tional fiberglass epoxy composite in Fig. 1, it has been shown that 
the angle of reflection of the reflected SV wave is smaller than that 
of the reflected P wave. Thus, the critical angle phenomenon does not 
occur for the unidirectional fiberglass epoxy composite under consid- 
eration. This is due to the fact that when the angle of incidence of 
an incident P wave reaches 90°, the angle of reflection of the 
reflected SV wave is still less than 90° 


4 Amplitude Ratios of Reflected Waves to Incident Wave 


It has been shown that when a P wave travelling in a semi- 
infinite transversely isotropic medium is incident on a plane bound- 
ary, a P wave and an SV wave will be reflected. The boundary condi- 
tions on the stresses of a P wave travelling in the plane x — 0 and 
incident on the plane boundary can be obtained from Eqs (4) and (7) 
as 


r yz 

r zz 


(I) 

(I) 


- 'yz (SV) 


+ r < sv > 

T r zz 


+ r yz 


+ r 


zz 


(P) - 0 
<p> - 0 


( 21 ) 
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where and r zz ^) represent the shear and normal stresses 

associated with the incident P wave; fy Z ^^ and r zz (^V) represent 
the shear and normal stresses associated with the reflected SV wave, 
and ty Z ^ p ) and r zz ^^ represents the shear and normal stresses 
associated with the reflected P wave. 


The shear stress, associated with a P wave of unit 

y i. 

amplitude travelling in the plane x - 0 and incident on the plane 
boundary at the origin, as shown in Fig. 2, can be obtained from 
Eqs . (3) and ( 6 ) as 


»yz (I) “ iu ( c 55 s z (I) 




( 22 ) 


where Sy(^) and S z ^) are the components of the slowness vector of the 
incident P wave along the y and z axes, respectively; Py(^) and P z ^ 1 ' 1 
are the components of the unit vector of particle displacement of the 
incident P wave along the y and z axes, respectively; C 55 is an elas- 
tic constant, and o> is the frequency. Similarly, the normal stress, 
r zz (^) , associated with the incident P wave of unit amplitude can be 
expressed, from Eqs. (3) and ( 6 ), as 


r zz (I ) - WC^SyWPyW + C 33 S Z C I >P Z (I>) (23) 

where C ^ 3 and C 33 are elastic constants. 
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on- 


The shear stress, r yz ( p ), and the normal stress, r zz ( p ), 
associated with the reflected P wave on the plane boundary at the 
gin are, from Eqs. (3) and (6), 


r y2 < P > - t U A< P ><C 55 S 2 < p >P y < p > + C 55 S y < p )p 2 ( p )) 

and 

r zz (?) - io,A(P)(C 13 S y (P)p y (P) + C 3 3S Z (P)P Z (P)) 


(24) 

(25) 


where S y ( p ) and S z ( p ) are the components of the slowness vector of the 
reflected P wave along the y and z axes, respectively; P y ( p ) and 
P z ^ p ) are the components of the unit vector of particle displacement 
of the reflected P wave along the y and z axes, respectively, and 
A< P ) is the amplitude of the reflected P wave. 


The shear stress, r yz ( sv ), and the normal stress, r zz ( sv ), 
associated with the reflected SV wave on the plane boundary at the 
origin, from Eqs. (3) and (6), 


r yz (SV) “ i<M/ SV >(C 5 5 S z ( SV )p y ( sv > + C 55 S y ( sv >P z ( sv >) (26) 

and 


r zz (SV) ” iwA( sv ) (C 13 S y ( sv )p y ( sv ) + C 33 S z < s V)p z ( sv >) (27) 
where S y ( sv ) and S z ( sv ) are the components of the slowness vector of 
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the reflected SV wave along the y and z axes, respectively; P y ( sv ) and 
P Z (SV) are t ^ e components of the unit vector of particle displacement 
of the reflected SV wave along the y and z axes, respectively; and 
a(SV) £ s t ^ e amplitude of the reflected SV wave. 

Upon substitution of Eqs. (22) through (27) into Eq. (21), the 
boundary conditions on the stresses for an incident P wave of unit 
amplitude travelling in the plane x - 0 can be rewritten as 


S z (I) Py ( 1 ) + Sy <I>P z <I> + A< p )(S z < P >P y < P > + S y < P >P z < P >) 

+ A(SV) ( s z (SV)p y (SV) + Sy (SV)p z (SV) ) = o 

(C 13 S y (I) P y (I) + c 33 s z ( I )p z ( 1 )) 

+ A( p )(C 13 S y ^ p )p y ( p ) + C 33 s z ( p )p z ( p )) (28) 

+ A(SV) ( c 13 s y (SV) Py (SV) 

+ c 33 s z ( sv )p z ( sv )) - o 

The components of a unit vector of particle displacement along 
the y and z axes of a P wave travelling in any plane containing the 
zonal axis of a transversely isotropic medium are given in [3], and 
can be expressed, when applied to the present case of the reflected P 
wave travelling in the plane x - 0, as [3] 
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( 29 ) 


P y (P) - [Hp - (C 33 -C 44 )S z ( p ) 2 ]/ 

( tH p -(C 3 3-C 44 )S z (P) 2 ] 2 

+ [ (C 13 +C 44 )S y ( P )S Z ( P ) ] ) ^ 

P z ( p) - (C 13+ C 44) s y ( p )s z ( p )/ 

( [Hp- (C 33 -C 44 ) S Z (P) 2 ] 2 

+ [ (C 13 +C 44 )S y ( P )S z ( P ) ] ) 

where Cjj, C33 and C 44 are elastic constants; and H p is defined as 

H p - ((C n -C 44 )S y < p > 2 + (C 33 -C 44 )S z ( p ) 2 

2 

+ {[(c n -c 44 )s y ( p > 

+ (C 33 -C 44 )s z ( p ) 2 ] 2 

- 4 S y ( p ) 2 S z ( p ) 2 [ (C n -C 44 ) (c 33 -c 44 ) 

- (Ci3+C 44 ) 2 ] 2 ) 1 /2 )/ 2. 

Similarly, the components of a unit vector of particle displacement 
along the y and z axes for the reflected SV waves travelling in the 
plane x - 0 can be expressed as 
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Py (SV) - [H SV - (C 33 -C 44 )S Z ( SV ) 2 ]/ 

{ [Hgv - (C 33 -C 44 )S z (SV) 2 ] 2 

+ [(C 13+ C 44 )S y (SV)s z (S v )] 2 } 1/2 (30) 

P z (SV) _ (C 13 +C 44 )S y < SV )s z < SV >/ 

{[H SV - (C 3 3-G44) S z (SV)2 ] 2 

+ [(C 13+ C 44 )S y (SV) Sz (SV) ) 2 ) 1 /2 

where H<;y is defined as 

H SV “ {(Cn-C 44 )S y < SV ) 2 + (C 33 -C 44 )S Z < SV ) 2 

- {[(Cn-C 44 )S y (SV) 2 + (C 33 -C 44 )S z (SV) 2 ] 2 

- 4S y < SV ) 2 S Z ( SV ) 2 [(C n -C 44 )(C 33 -C 44 ) 

2 2 1/2 

- ( C 13 +C 44) ] ) 5/2 • 

The amplitude ratios of the reflected P and SV waves, A^P) and 
A< sv >, for a P wave of unit amplitude obliquely incident on the plane 
boundary can be determined from Eq. (28). For a given incident P 
wave, the values of the components of the slowness vector along the y 
and z axes, and s z ^', and the values of the components of the 

unit vector of particle displacement along the y and z axes, P y ^5 and 
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P z ^), are defined as part of the the specification of the incident P 
wave The values of the components of the slowness vector of the 
reflected P wave along the y and z axes, Sy( p ) and are 

determined from Eqs . (16) and (17), respectively. The values of the 
components of the unit vector of particle displacement along the y and 
z axes for the reflected P wave, Py'^) and P z ^) , are obtained by 
substituting Sy^^ and into Eq. (29). The values of the 

coefficients for in Eq. (28) are thus obtained Since the y-com- 

ponent of the slowness vector of the reflected SV wave, , is 

equal to that of the incident P wave, Eq. (16), the value of the 
z-component of the reflected SV wave, S Z ^ V ), is obtained by substi- 
tuting Sy(^) into the slowness surface for the SV wave, Eq. (15). On 
substitution of the values of the S z ^^^ and Sy(^V) i n £ 0 Eq. (30) , the 
values of the unit vector of particle displacement along the y and z 
axes for the reflected SV wave, Py(SV) anc j p^SV), are determined 
The values of the coefficients for A^ v ) in Eq. (28) are thus 
obtained. The values of the amplitudes of the reflected P and SV 
waves , A<*> and A<SV>, are then obtained by solving Eq. (28) with the 
thus determined values of the coefficients for A^) and A^^ in Eq 
(28) for a given incident P wave. By varying the angle of incidence 
of the incident P wave of unit amplitude and by repeating the proce- 
dures described above, the amplitude ratios of the reflected P and SV 
waves to the incident P wave are obtained as functions of the angle of 
incidence. For an incident P wave in the unidirectional fiberglass 
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epoxy composite shown in Fig. 1, the amplitude ratios of the reflected 
P and SV waves to the incident P wave versus the angle of incidence 
are shown in Fig. 5. 

5 Balance in Energy Flux Normal to Plane Boundary 

The balance in energy flux normal to the plane boundary z = 0, 
as shown in Fig. 2, must be satisfied and is expressed as [2] 

F z d) + F Z (P) + F Z < SV > - 0 (31) 

where F z ^) , F z ^^ and F Z ^V) are the z-components of the energy 
fluxes of the incident P wave, the reflected P wave and the reflected 
SV wave, respectively. 

The z- component of the energy flux of an incident P wave of unit 
amplitude travelling in the plane x - 0 is [4] 

F Z <I> - „2 (044 p y <I> 2 s z «) ♦ C 13 P y (I)p 2 (I). y (I) 

+ C 44 P y <I>P z <»S y «> + C3 3 P z (I) 2 S z (I), (32) 

where Py^ 1 ) and P z ^ are the components of the unit vector of 
particle displacement of the incident P wave along the y and z axes, 
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respectively; S y (*) and S 2 ^) are the components of the slowness 
vector of the incident P wave along the y and z axes, respectively, u> 
is the radian frequency, and C 33 , C 33 and C 44 are elastic constants 

Similarly, the z-components of the energy fluxes of the 
reflected P wave and the reflected SV wave are [4] 

F Z < P > - A<P > 2 U 2 (C 44 P y < p > 2 S Z < P > 4- C 13 P y < p >P z < p >S y < P > 

+ C 44 p y < p )p z < p >s y < p > + C 33 P z < P > 2 S z < p >) (33) 

and 

F Z (SV) _ A ( SV > 2 W 2 (C 44 P y (SV) 2 s z (SV) + c 13 P y (SV)p z (SV) Sy (SV) 

+ C 44 P y ( SV )p z ( SV >S y ( SV ) + C 33 P z ( SV ) 2 S z < SV >) (34) 

Since the frequency term w 2 is common to Eqs . (32), (33) and 
(34), the balance in energy flux normal to the plane boundary, 

Eq (31) , is not affected by assuming the value of the frequency o> to 
be equal to unity. Accordingly, subsequent calculations of the values 
of the z-components of the energy fluxes of the incident P wave and 
the reflected P and SV waves, F z ^ , F Z ^ P ^ and F Z ^ SV \ are done by 
assuming the radian frequency w to be equal to one. 
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The values of the z- components of the energy fluxes in 
Eqs. (32), (33) and (34) are obtained similarly to the calculations 
of the amplitude ratios of the reflected P and SV waves to the inci- 
dent P wave. For a given incident P wave of unit amplitude, the 
y-components of the slowness vectors and of the unit vectors of par- 
ticle displacement, SyX*) , Sy( p ) , Sy( sv ) , Py(*), Py( p ) and Py( sv ) , and 
the z-components of the slowness vectors and of the unit vectors of 
particle displacement, S z ^^, S z ^^, S Z ^ V ), P z ^^, P Z ^ P ^ and P z ( sv ) 
of the incident P wave and of the reflected P and SV waves are deter- 
mined first. Then, combining the amplitudes of the reflected P and SV 
waves, and , with the values of Sy^^ , Sy(^) , Sy^ v ) , s z ^\ 



values of the z-components of the energy fluxes of the incident P wave 
and the reflected P and SV waves, F z ^ ^ , F z ^ p ^ and F z ^ sv \ are thus 
obtained from Eqs. (32), (33) and (34). For a P wave of unit ampli- 
tude travelling in the plane x - 0 in the unidirectional fiberglass 
epoxy composite shown in Fig. 1 incident on a plane boundary, the 
z-components of the energy fluxes of the reflected P and SV waves and 
the energy flux of the incident P wave are shown in Fig. 6 with the 
value of the frequency w in Eqs. (32), (33) and (34) equal to one The 
balance in energy flux normal to the plane boundary, Eq. (31), is also 
checked and is shown in Fig. 6. 
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ACOUSTO-ULTRASONIC NON-CONTACT INPUT-OUTPUT CHARACTERIZATION OF UNIDI- 


RECTIONAL FIBERGLASS EPOXY COMPOSITE PLATE 


It has been shown [1] that the unidirectional fiber composite 
shown in Fig. 1 may be modelled as a homogeneous transversely iso- 
tropic continuum For the axes shown in Fig. 1, the isotropic plane 
of its equivalent continuum lie in the midplane of the plate [1], A 
cartesian coordinate system (x, y, z) is chosen so that the x-y plane 
is the isotropic plane; as a result, the upper and the lower surface 
are at z - h/2 and z - -h/2, respectively, where h is the plate thick- 
ness The properties of the equivalent continuum model of the unidi- 
rectional fiberglass epoxy composite plate to be considered are [1] 

h - 0.1 m 

C n - 10.581 x 10 9 N/m 2 

C 13 - 4.679 x 10 9 N/m 2 

C 33 - 40.741 x 10 9 N/m 2 (35) 

C 44 - 4.422 x 10 9 N/m 2 

C 66 - 3.243 x 10 9 N/m 2 

p - 1850 kg/nr* 

A transmitting and a receiving transducer are located on the 
same face of the unidirectional fiberglass epoxy composite plate spe- 
cimen without direct contact, as shown in Fig. 7. The unidirectional 
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fiberglass epoxy composite plate specimen shown in Fig. 7 is consid- 
ered as a plate of thickness h and of infinite planar (x-y) extent 
The input electrical voltage to the transmitting transducer is V^(t) 
and the output electrical voltage from the receiving transducer is 
V Q (t) where t represents time. The transmitting transducer converts 
an input electrical voltage into a stress, whereas the receiving 
transducer converts a displacement associated with stress waves trav- 
elling in the plate into an output voltage. In the following analy- 
sis, only the P waves are traced. The P waves which are generated by 
the transmitting transducer located above point 0 experience multiple 
reflections at each face of the plate, and then reach the receiving 
transducer located above point M, as shown in Fig 8. Since the iso- 
tropic plane lies in the midplane and is parallel to both the top and 
the bottom faces where the multiple reflections occur, the angle of 
reflection of the reflected P wave is equal to the angle of incidence 
of an incident P wave for each reflection at each face of the plate 
Accordingly, the P wave travelling from the input 0 to the output M 
may be considered as a wave propagating in a semi- inf inite trans- 
versely isotropic medium, and travelling to point M' as if there were 
no bottom face, as shown in Fig. 8. 

1 • Delay Time and Phase Velocity 

Let the input 0 and the output M lie in the y-z plane. Assume 
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the number of reflections at the bottom face experienced by the P 
wave in travelling from the input 0 to the output M is n, as shown in 
Fig. 8 . With respect to the z axis, the angle of incidence of the P 
wave at each face of the plate is 9, and the total distance travelled 
by the wave is E^. From the geometry in Fig. 8 , 

9 - tan' 1 (i/2nh) (36) 

where SL is the separation distance between the input 0 and the output 
M, and 

E^ - i/sin 9 (37) 

The time delay t n for the wave to reach the receiving transducer is 

t n “ Rn/ C l(*) (38) 

where C]_(0) is the directionally dependent phase velocity of the P 
wave. The phase velocity of a P wave in the unidirectional fiber- 
glass epoxy composite shown in Fig. 1 is given as [1] 

Ci(9) - [(C 44 +C 1;L sin 2 5 + C 33 cos 2 9 + J^/lp] 1 / 2 (39) 

where 

E - [(C U -C 44 ) sin 2 0 + (C 44 -C 3 3 )cos 2 5] 2 
+ ^(C^ 3 +C 44 ) 2 sin 2 0 cos 2 0 ; 
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c ll> c 13> c 33 > c 44 and P are given by Eq. (35). 


The delay time is then computed when the number of reflections n 
at the bottom face of the plate is equal to 10, 100, 300 or 500 The 
numerical results are shown in Fig. 9 where the delay time t n is 
plotted as the ordinate, and the dimensionless separation 2 / h is 
plotted as the abscissa, for values of zero to 300. The phase velo- 
city as a function of the angle of incidence 6 is also shown in 

Fig 10. 

2 . Displacements Detected by Receivine Transducer 

The displacements detected by the non-contact receiving trans- 
ducer above point M, radiated by the non-contact transmitting trans- 
ducer, are assumed to be equivalent to the displacements at point M' 
associated with the P wave propagating in a semi- infinite transversely 
isotropic medium as if there were no bottom boundary (except for the 
cumulative effect of the reflection coefficient), as shown in Fig. 8 
The displacements at point M' are approximated by the far- field asymp- 
totic solution for large of an infinite transversely isotropic 
medium subjected to a harmonic point load. 

Consider an infinite transversely isotropic medium in which the 
z axis of a rectangular cartesian system 0(x, y, z) is the zonal axis 
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of the medium and the x-y plane coincides with the isotropic plane, as 
shown in Fig. 11. The equations of motion including the body forces 
are [5] 

r xx, x + T xy ,y + r xz,z + P X “ pu.tt (40) 

r xy , x + r yy,y + r yz,z + pY - pv.tt (41) 

T xz,x + r yz ,y + r zz,z + P z “ P w * t,: ( 42 > 

where r rs (r, s-x, y and z) are the normal (r — s) and shear (r ^ s) 
stresses with respect to the chosen coordinate system 0(x, y, z) ; u, 
v and w are displacement components of a point in the medium along the 
x, y and z axes, respectively; X, Y and Z are the components of the 
body force along the x, y and z axes respectively; p is the density, t 
is time; and denotes partial differentiation with respect to the 
variable which follows. 


Combining Eqs. (l)and (42) gives 


r.tt- <- 


c 13 +c 44 


-44 


-44 


P 

c 33 


} a 4- r 4- — r 

' zz ^ p p i, yy 


^ ' ZZ X > ; 


where T and A are given by [ 5 ] . 


(43) 
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r - w,z; A - u,x + v,y . 


By differentiating Eq. (40) with respect to x and Eq. (41) with 
respect to y, we find upon addition of the resulting equations and 
using the appropriate stress-strain relations, Eq (1), that 


A, tt 


C13+C 


44 C44 

) (r,xx + T.yy) +— — A,zz 


'11 


(A, xx + A.yy) + X,x + Y,y 


(44) 


For a harmonic point load at the origin, the body forces may be 
taken of the form [5] 

X - X 0 5(x)5(y)5(z)e* i “ t 

Y - Y 0 S(x)S(y)5(z)e* lwt (45) 

Z - Z 0 S(x)6(y)S(z)e- iwt 

where 6(r) (r - x, y and z) is the Dirac delta function, and X Q , Y 0 
and Z Q are the magnitudes of the respective point body forces. 
Express T and A as threefold Fourier integrals [5] as follows 
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CO 


T(x,y , z , t) 


r ( s x , Sy , s z , t) 


-co 


exp{iw(S x x+S y y+S z z-t) }dS x dS y dS z 


( 46 ) 


A(x,y,z,t) 



A(S X , S y , S z , t) 


exp(iw(S x x+S y y+S z z-t) )dS x dS y dS z 


where 


r(S x ,S y ,S z ,t) - 1/8jt' 


r(x,y,z,t) 


/ W a 

-co 


exp{iw(-S x x-S y y-S z z+ t))dxdydz. 


A(S x ,S y ,S z ,t) - l/8w 3 



-CO 


A(x,y,z,t) 


exp{iw(-S x x-S y y-S z z+ t))dxdydz. 


Similarly, X,x, Y,y and Z,z can be expressed as threefold 
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Fourier integrals [5] 


co 


X,x - iwS x 


X exp{iw(S x x + Syy + S z z - t))dS x dSydS z 


-CO 


(47) 


Y.y - iws. 


Y exp{iw(S x x + 


s y y 


S z z 


t) JdS x dSydS z (48) 


Z,z - iwS. 


Z exp{iw(S x x + Syy 


S z z 


- t))dS x dS y dS z (49) 


-00 


where 


X 


1/8* 3 



X 0 S(x)S(y)S(z)e' iwt 


exp{iu(-S x x - Syy - S z z + t)dS x dSydS z 

- X 0 /8jt 3 
and, similarly, 



Substitution of Eqs. (45) through (49) into Eqs. (43) and (44) 
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gives 


C44 Jp (S x 2 +S 2 ) + C33/P S z 2 -1 

— - iL — — 5 (iS x X 0 + iS v y o ) + 

O— C O N A.U yu" 


8jr J uH (S x , S yf S z ) 


( c 13+ c 44 )/p (S x 2 +S y 2 ) 

8ir 3 cH (S x , S y , S z ) 1 Z ° 


( 50 ) 


-(C 13 +C 44 )/p S z 


-3 — - — — (iS x X 0 + iS y Y 0 ) 

8* 3 wH (S x , S y , S z ) X 0 y 0 


C 11 /P (S x 2 + S y 2 ) + C33/P S z 2 - 1 _ 


8 tt j uH (S x , S y , S z ) 


is z z o 


(51) 


where 


H(S x ,S y ,S z ) 



(S x 2 +S y 2 ) 



c 44 


(S x 2 +S y 2 ) + 





c 44 +c 13 s 2 „ 2 

z ) i> z 



(52) 


In fact, H(S X , S y , S z ) - 0 represents two sheets of the slowness 
surface, one for a P wave and one for an SV wave [5] As a result of 
Eqs . (50) and (51), T and A in Eq. (46) can be written as 
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r (x,y , z , t) - 


-( C 13 + C 44 )/P S z 2 

— (iS x X 0 +iS„Y 0 ) 

8*\>H(S X ,S S z ) X 0 y 0 


4 cn /p (s x 2 +s v 2 ) + c 33 / P s z 2 - l 

~ 8A>H(S x ,S y ,S z ) 1 Z ° 

exp { iu(S x x+S y y+S z z-t) )dS x dS y dS z 


( 53 ) 


A(x,y,z,t) 


CvJe + Z Ji/l S " 2 ~ 1 ( iS x X o +lS y Y 0 ) 

8* 3 u>H(S x ,S v ,S 7 ) x 0 y 0 


■»x» u y » ,J z> 

, -(Cj3 + C44)/P (S x 2 +S y 2 ) ' 

8fl’ 3 wH(S x , S y , S z ) 1 z 0 


exp {iw(S x x + S y y + S z z - t)}dS x dS y dS z 


(54) 


The asymptotic solution at a large distance from the point load is 
obtained by applying the theory of residues, the method of stationary 
phase, and the radiation condition [ 5 ] as 


r(x,y , z , t) - 


A n“ 

2?rR 


•(C13+CAA.) *2, * * 

S z Z (iS x X Q +iS Y 0 ) 


+ (Sx 2+ s* 2 ) + -^- 3 S * 2 - l)iS* z 0 


exp{iw(S*x + S*y + S*z - t) ) 


A(x,y,z,t) 


f 

2?rR \ 


c 44 ,„*2. „*2, 
P 


(55) 


Hr-( V+ s/) + 


S* 2 - l](iS* X 0 +iS y Y 0 ) 


-(C13+C44) *0 *2. * 

+ — — =-*— (S./+ S z ) iS_ z, 

p X y'z 


•} 


exp{iw(S x x + S y y + S z z - t) ) 


(56) 
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where R is the distance from the origin 0 where the point loads are 

applied to the location of interest Q in the medium, as shown in 
* * * 

Fig. 11; (S x , Sy, S z ) are points on the slowness surface for a P wave 
where the normal is parallel to the OQ direction; A n is the amplitude 
coefficient and is given by 



where 2 is the sum with respect to cyclic permutation of S x , Sy and 

* * * 

S z , and is evaluated at points (S x , Sy, S z ) on the slowness for a P 
wave where the normal is parallel to the OQ direction. 


The displacement components along the x, y and z axes, u, v and 
w, due to a P wave can be obtained by direct integration of the defi- 
nitions of T and A in Eq. (43), and are given as [6] 

s x 

u - A 

iw(S x 2 +S y 2 ) 

s y 

v A (58) 

iO)(S X 2 +Sy 2 ) 

1 

w r 

iu>s z 
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Substitution of Eqs. (55) and (56) into Eq. (58) gives the 
asymptotic solutions at a large distance of the displacement compo- 
nents along the x, y and z axes, u, v and w, as follows: 


•i*n S* 


2ttR(S* 2 + S* 2 ) 


C44 *2 *2. C33 *2 * 

Hr- (s x + s y > + rr s z -lias* x o + 


iSy Y q ) 


-(C 13 +C 44 ) * 2 . „*2 


(s x + s v ) is z Z Q 


^ 

exp{iw(S x x + Syy + S 2 z - t) ) 


(59) 


v - 


■Un 4 


2jtR(S* 2 + S* 2 ) 


r C 44 /n *2 . e*2 s , ^33 „*2 

P 


p£- (s "‘ + S^) +-^ S z *- l](iS x X 0 + iS* Y 0 ) 


+ :< C 13^44) (s; »2 + S *2 )1S * , 




ic ic ic 

exp{iw(S x x + Syy + S z z - t) } 


(60) 


w - 


’n 


2rrR 


-( c 13 + C 44 ) * * * 

“ ^ - S z <s x X 0 + Sy Yq) 


♦ A* (.?♦»?> i)z t 


"ft "Jc 'fc 

exp{iw(S x x + S y y + S z z - t) ) 


(61) 
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The slowness surface for a P wave travelling in the unidirec- 
tional fiberglass epoxy composite shown in Fig. 1 is an oblate sphe- 
roid See Fig. 3. Thus, if (x, y, z) are the coordinates of a given 

point in the medium, there will be only one P wave front passing 

* * * 

through it, corresponding to a point (S x , Sy, S z ) on the slowness 
surface where the normal is parallel to the given direction. There- 
fore, it is concluded that for the displacement components along the 
x, y and z axes, u, v and w, measured by tracing the P wave in the 

unidirectional fiberglass epoxy composite plate specimen at point M m 

•k * 

Fig 8, there is only one point (S x , Sy, S z ) on the slowness surface 
where the normal is parallel to a given direction OM' [5]. 


3 . Directivity Function 

The shear stress r y Z and the normal stress r zz associated with 
the P wave reaching the point M' in Fig. 8 are used to study their 
associated directivity functions. The asymptotic shear stress r y Z and 

the asymptotic normal stress r zz are obtained by, substituting 

■* 

Eqs (55) through (58) into Eq. (1) and then by setting S x equal to 
zero as 
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zz 


A n w f -( c 13 +c 44) * *2 v 

2jrR C 33 p LS y s z Y o 


, / C H o*2 , C 33 *2 .. . * . r c 44 *2 

+ ( p s y + p s z ' is z z oJ + c 13 t (~ “ ^y 


, c 33 „*2 1S „ s . *(Ci3+C 44 ) , o * 2o * „ ,, 

+ p ^z*l) L ^y V + p lS y s z z o]} 


exp{iw(Syy + S z z - t) } 


( 62 ) 


c 55 A n“ f -(C 13 +C 44 ) ,„*2„* „ 

r y z " { — p — 1 y Sz Y ° 


'll *2 , ^33 „*2 


+ ( S v + 

p y p 


S Z - 1) iSy Z Q + (- 


J 44 *2 

~p~ b y 


^ s : 2 - x, 14 v ♦ - <Ci3 ; C44) is* s: 2 z o) 


exp{ioj(Syy + S z z - t) ) 


(63) 


The directivity functions associated with the normal stress r zz 
in Eq. (62) and the shear stress Ty Z in Eq. (63) will be evaluated 
Due to the axial symmetry with respect to the zonal axis, the z axis, 
of the transversely isotropic medium, the values of the directivity 
functions obtained for the case S x - 0 hold for all values of S x 


39 



Consider the case of the point load acting along the y direction 
only; that is, Y 0 ** 0 but X Q - Z Q - 0 in Eq. (45). 

Y 

The directivity function Dy Z of the shear stress r y Z associated 
with a P wave whose slowness vector is confined to the plane x - 0 is 
obtained from the amplitude of the shear stress in Eq (63) by setting 
Y 0 - R - 1; 


* 

Y A n wC 44 S z f -C u ^*2 c 33 „*2 , „ 

)yz — (-y- S y +— b z - 1) 


(64) 


where A n is given in Eq. (57) and is evaluated at the point (0, Sy, 

S z ) on the slowness surface for a P wave, Eq (14); (0, Sy, S z ) is the 

point on the slowness surface where the normal is parallel to the 

given direction and u is radian frequency. Similarly, the directivity 
Y 

function D zz of the normal r zz is, from Eq. (62), for Y 0 - R - 1, 



■X n o>S y 

2 * 


c 13 c 44 „*2 c 33 c 44 

“ — — — 



(65) 


Next, consider the case of the point load acting along the z 
direction only; that is, Z Q * 0, but X 0 - Y 0 - 0 in Eq. (45). 


Z 

The directivity function Dy Z of the shear stress Ty Z associated 
with a P wave whose slowness vector is confined to the plane x - 0 is 
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obtained from the amplitude of the shear stress in Eq. (63) by setting 


R - 1; 


* 

z A n wC 44 S y Cn o * 2 C 33‘C 44 *C 13 „*2 

yz ” „ b y + b z ‘ 

J 2n P J p 


( 66 ) 


Z 

Similarly, the directivity function D zz of the normal stress r zz is, 
from Eq. (62), for Z 0 - R - 1; 


★ 2 2 
A n wS z , c ll c 33- c 13 c 44- c 13 „*2 . c 33 ^*2 ^ N /r , x 

D zz ( s y + s z * c 33^ ( 67 > 

2n p J P 

The polar diagrams for the directivity functions of the shear 

Y Y Z Z 

stress Ty Z and the normal stress r zz (Dy Z , D zz , Dy Z and D zz given in 
Eqs. (64) through (67)) associated with the P waves propagating in the 
unidirectional fiberglass epoxy composite shown in Fig. 1 are obtained 
by substituting Eq. (35) into Eqs. (64) through (67) at frequencies of 
0 75, 1.50 and 2.25 MHz. 

Numerical results are shown in Figs. 12 through 23 where the 

angle of incidence 6 given in Eq. (36) is used to determine the 

* * 

direction for which the value of a point (0, Sy, S z ) on the slowness 
surface where the normal is parallel to the given direction is thus 
obtained. 


41 



y. Assumptions on the Transducers 

The non- contact transmitting transducer in Fig. 7 is assumed to 
transform an electrical voltage into a uniform stress; however, the 
non-contact receiving transducer in Fig. 7 transforms a displacement 
into an electrical voltage. The approach below is similar to that 
given in [7] . 

Referring to Fig. 7, if an input voltage of amplitude V and 
frequency w is applied according to 


V^t) - Ve' iwt 


( 68 ) 


the stress a that is introduced into the specimen plate by the 
non- contact transmitting transducer is 

c(t) - F 1 (u)Ve- i ( wt + ^1> (69) 

where F^(w) is the transduction ratio for the non-contact transmitting 
transducer in transforming a voltage to a stress and <j > i is a phase 
angle. In Eqs. (68) and (69), the harmonic character of the signals 
is expressed in the complex notation where i- J- 1 and only the real 
parts of these and subsequent equations should be considered. Thus, 
the amplitude T of the applied force is defined as 
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T - 


(70) 


Similarly, if a stress wave producing a displacement d of ampli- 
tude D and frequency w that is detected by the non-contact receiving 
transducer is defined as 

d(t) - De- iwt (71) 

the output voltage from the non-contact receiving transducer, see 
Fig. 7, is 

-i(wt + <f> 2 ) 

V Q (t) - F2(w)De (72) 

where F 2 (w) is the transduction ratio for the non-contact receiving 
transducer in transforming a displacement to a voltage, and <j> 2 is a 
phase angle Thus, the amplitude V' of the output electrical voltage 
is 


V' - F 2 (w)D 


(73) 


The characteristics of F^(w) and F 2 (cj) 
the dimensions of the product F^(w)F 2 (w) are 


are unknown except that 

[kg/m^'sec^] . 
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tions In a Plate 


Since the P wave traced in the unidirectional fiberglass epoxy 
composite plate specimen shown in Fig. 8 is travelling in the y-z 
plane, it follows from Eqs. (59), (60) and (61) that only the dis- 
placement components along the y and z axes, v and w, are detectable 
at the point M' Consider first the point load acting along the y 
direction only; that is, Y Q ** 0, but X 0 - Z Q - 0 in Eq. (45). 


Y 

The amplitude of the y-component displacement D v evaluated at 
the point M' can be obtained from Eqs. (37) and (60) as 


* * 

Y fl(Sy>Sz) 

D v - Y 0 (74) 

*n 


* ★ *2 *2 
where f^ (Sy,S z ) - (A n /2»r) | (C 44 /P) S y + (C 3 3 /V)S Z -l|. Similarly, the 

Y 

amplitude of the z-component displacement D w evaluated at the point M' 
can be obtained from Eqs. (37) and (61) as 


* * 

Y 

D w Y 0 (75) 

*n 


* * * * 

where f 2 (S y ,S z ) - (A n /2*) [ (^ 3 + 044 )/?] S y S z . 
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Next, consider the point load acting along the z-direction only, 
that is, Zq r* 0 , but X Q - Y Q - 0 in Eq. (45). 

Z 

The amplitude of the y-component displacement D v evaluated at 
the point M' can be obtained from Eqs . (37) and (60) as 

* * 

Z f 2( s y> s z) 

D v L z o (76) 

*n 

evaluated at the point M' can be obtained from Eqs. (37) and (61) as 
* * 

Z f 3( s y> s z) 

D w L z o (77) 

Rn 

where f 3 (S*, S*) - (A n /2»r) | (C 11 /p)S* 2 +(C 33 /p)S* 2 -l | . 

Finally, consider the point load acting along the x direction 
only; that is, X Q * 0, but Y Q - Z Q - 0 in Eq. (45). It follows from 
Eqs. (60) and (61) that the amplitude of the y-component displacement 
and the amplitude of the z-component displacement vanish at the point 
M' . This is due to the fact that the P wave is travelling in the y-z 
plane . 

According to Eqs. (74) through (77), the amplitude of a displa- 
cement component at the point M' , denoted as D^» , can be expressed in 
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the following form as 


* * 
fi(Sy,S z ) 

D m , T (78) 

*n 

where T is the applied point load and is equivalent to either Y 0 or 

★ ★ 

Z Q in Eqs. (74) through (77); and f£(Sy,S z ) (i-1,2 or 3) is determined 
from one of the Eqs. (74) through (77), depending on which displace- 
ment component is measured and along which direction the point load is 
applied. 

However, with the bottom boundary present, the wave is reflected 
a total of (2n-l) times, as shown in Fig. 8. Thus, the amplitude of 
displacement at the point M is and is expressed as 

2n-l 

d M “ Qpp d M' ( 79 > 

where Qpp is the amplitude ratio of the reflected P wave to the 
incident P wave, as shown in Fig. 5, and is a function of the angle of 
incidence 9 defined in Eq. (36). 

Here, we have ignored the effects of mode conversion. The ampli- 
tude of the output voltge from the non-contact receiving transducer is 
V' and can be obtained by substituting Eqs. (78) and (79) into Eq 
(73) as 
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V' 


(80) 


2n-l * * 

f'2( CJ ) Qpp f i(Sy , S Z )T 

Rn 

Substitution of Eq. (70) into Eq. (80) gives 

2n-l * * 

F i(w) F 2 («) Qpp fi(S y , S Z )V 
V' — - (81) 

Rn 

Introducing the P wave attenuation constant a of the unidirectional 
fiberglass epoxy composite and a possible electrical signal amplifica- 
tion factor K, Eq. (81) can be written as 

V' - KF 1 (w)F 2 (w)Qpp" 1 f i (Sy , S*)V — (82) 

Eq. (82) gives the output voltage amplitude from the non-contact 
receiving transducer due to an *nput voltage amplitude V at the non- 
contact transmitting transducer when the P wave path has undergone n 
reflections from the bottom face of the unidirectional fiberglass 
epoxy composite plate specimen, as shown in Fig. 8. 
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DISCUSSION AND CONCLUSION 


In the acousto-ultrasonic input-output characterization of the 
unidirectional fiberglass epoxy composite plate specimen, the angle of 
reflection of the reflected P wave is equal to the angle of incidence 
of the incident P wave for each reflection at either the top or the 
bottom face of the plate. This is due to the fact that the isotropic 
plane is parallel to both faces of the plate However, if the paral- 
lelism between the isotropic plane and the plane boundaries where 
reflection occurs does not exist, the angle of reflection is not equal 
to the angle of incidence. In this case, the use of a semi- inf inite 
transversely isotropic medium, neglecting the existence of the bottom 
face of the plate specimen except for the cumulative reflections, to 
compute the delay time, the displacement and the directivity functions 
becomes inappropriate. 

It is observed from Fig. 9 that an increase in the number of 
reflections n results in a minor increase in the delay time t n at 
each value of the dimensionless separation i/h. For a given plate 
thickness h and a given the separation distance i between the trans- 
mitting transducer and the receiving transducer, an increase in the 
number of reflections n results in a decrease in the angle of inci- 
dence 9 , as given by Eq. (36), thereby increasing the travelling 
distance from Eq. (37). The phase velocity of a P wave in 
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the unidirectional fiberglass epoxy composite increases with decreas- 
ing angle of as shown in Fig. 10. For the time delay, the increase 
in the travelling distance is partly offset by the increase in the 
phase velocity with an increse in the number of reflections n 
This accounts for the minor increase in the delay time given by Eq 
(38). However, in an isotropic medium, the phase velocity of a stress 
wave is directionally independent. As a result, the increase in the 
dealy time t n caused by an increase in the number of reflections n is 
solely attributed to the increase in the travelling distance in an 

isotrpic plate. Therefore, it is concluded that the directional 
dependence of the phase velocity of a stress wave travelling in a 
transversely isotropic medium has a significant effect on the delay 
time when conducting acousto -ultrasonic testing. 

Now, consider the case of a given number of reflections n. An 
increase in the separation distance £ results in an increase in the 
travelling distance R n , as shown in Fig. 8, and an increase in the 
angle of incidence 0 given by Eq. (36). An increasing angle of inci- 
dence 0 results in a decreasing phase velocity C^, as shown in 
Fig. 10. Consequently, the delay time t n increases sharply as a 
result of the increasing separation distance £ for a given number of 
reflections n, as shown in Fig. 10. However, as the separation dis- 
tance i approaches infinity, the phase velocity approaches a limit. 
Thus, (for i/h-wo) the increase in the delay time t n is due solely to 
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the increase in the travelling distance which is similar to the 

isotropic medium case. 

This theoretical investigation provides a step forward in the 
quantitative understanding of acousto-ultrasonic nondestructive evalu- 
ation (NDE) parameters such as stress wave factor (SWF) in trans- 
versely iostropic media. It also provides the potential for assisting 
in the development of more efficient and more revealing NDE schems 
utilizing wave propagation. 
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Fig. 2 Coordinate system (x,y,z) in analysis of single reflection 
problem at stress-free plane boundary of semi-infinte 
transversely isotropic medium; z = 0 is plane boundary where 
single reflection occurs, z < 0 is free space, and x = 0 is 
plane of incidence. 
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Sz( s/km) 


Fig. 3 



Component of slowness vector in y direction, Sy(s/km) 


Two sheets of slowness surfaces for P and SV waves in unidirectional fiberglass epoxy 
composite for positive y-z quadrant. 
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Fig. 4 Angle of reflection of reflected SV wave versus angle of incidence of 
incident P wave in unidirectional fiberglass epoxy composite. 
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Fig. 6 Energy fluxes normal to plane boundary, assuming radian frequency equal to unity, 
for incident P wave m unidirectional fiberglass epoxy composite. 
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Fig. 7 Schematic of non-contact acousto-ultrasonic test configuration 
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Fig. 9 Delay time versus dimensionless separation ( £/h) with number of reflections n, 
from bottom face of fiberglass epoxy composite plate specimen, as parameter. 
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Fig. 12 


Polar diagram for directivity function D yz of shear stress T yz at transmitting 
origin associated with P wave in unidirectional fiberglass epoxy composite at 
frequency of 0.75 MHz for positive y-z quadrant due to applied point load acting 
along y-direction. 
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Fig. 15 Polar diagram for 
origin associated 
frequency of 0.75 
along y-direction 
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Fig. 16 Polar diagram for directivity function D zz of normal stress T zz at transmitting 
origin associated with P wave in unidirectional fiberglass epoxy composite at 
frequency of 1.50 MHz for positive y-z quadrant due to applied point load acting 
along y-direction. 
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Polar diagram for directivity function D yz of shear stress t at transmitting 
origin associated with P wave in unidirectional fiberglass epoxy composite at 
frequency of 1.50 MHz for positive y-z quadrant due to applied point load acting 
along z-direction. 
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Fig. 20 Polar diagram for directivity function D yz of shear stress T yz at transmitting 
origin associated with P wave m unidirectional fiberglass epoxy composite at 
frequency of 2.25 MHz for positive y-z quadrant due to applied point load acting 
along z-direction. 
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Fig. 23 Polar diagram for directivity function D Zz of normal stress T zz at transmitting 
origin associated with P wave in unidirectional fiberglass epoxy composite at 
frequency of 2.25 MHz for positive y-z quadrant due to applied point load acting 
along z-direction. 
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